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Abstract 

Some results in random matrices are generalized to supermatrices, in particular superma- 
trix integration is reduced to an integration over the eigenvalues and the resulting volume 
element is shown to be equivalent to a one dimensional Coulomb gas of both positive and 
negative charges.lt is shown that, for polynomial potentials, after removing the instability 
due to the annihilation of opposite charges, supermatrix models are indistinguishable from 
ordinary matrix models, in agreement with a recent result by Alvarez-Gaume and Manes. 
It is pointed out however that this may not be true for more general potentials such as for 
instance the supersymmetric generalization of the Penner model. 
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One would naively expect that a model based on supermatrices is a generalization of 
matrix models [1] suitable to describe a discretized two dimensional supergravity coupled 
with superconformal matter. In fact it contains both fermionic and bosonic degrees of 
freedom related by a local ( in a zero dimensional sense) super symmetry. 

However in matrix models the matrix elements are sort of preons whose interpretation 
in terms of two dimensional space time is far from being obvious. So the supersymmetry in 
supermatrix models may easily not have anything to do with space-time supersymmetry. 

We show in this note that this is indeed the case, and that, in supermatrix models 
with an action S polynomial in the supermatrix A, the fermionic degrees of freedom simply 
cancel an equal number of bosonic degrees of freedom, so that in the end one is left with an 
ordinary matrix model with a reduced (if any) number of entries. This confirms the results 
obtained in other way in ref. [2]*. We point out however that with non polynomial actions 
such cancellation between bosonic and fermionic degrees of freedom does not necessarily 
occur, leading possibly to models that are not equivalent ordinary matrix models. 

We show as well that, in the same way as ordinary matrix models are equivalent to a 
Coulomb gas of equal charges, whose positions are given by the eigenvalues of the matrix, 
supersymmetric matrix models can be described in terms of a Coulomb gas of both positive 
and negative charges whose positions correspond to the eigenvalues of the supermatrix. The 
sign of the charge depends upon which bosonic submatrix the corresponding eigenvalue 
belongs to. The annihilation of a couple of opposite charges corresponds in this picture to 
the cancellation of bosonic and fermionic degrees of freedom. 

Before considering the relatively more complicated case of matrix models let us study 
the supersymmetric generalization of the vector models analyzed in [3], where supersymme- 
try leads to a rather straightforward cancellation between bosons and fermions. Consider 
a vector v consisting of n bosonic and m fermionic components: 

/xi \ 



and a partition function given by 

Z = / We -« (!) 
where v 2 is the OSp(n, m) invariant scalar product: 

n m 
i=l a=l 

* Most results contained in this paper, and in particular eq.(13) have been obtained 
prior to our knowledge of ref. [2] 
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and C af3 is a symplectic matrix. 

One can compute Z by gauge fixing the OSp(n, m) symmetry. With the gauge fixing 
conditions 

X 2 = X 3 = ■ ■ ■ = X n = = 1p 2 = ■ ■ ■ = = 

the Faddeev-Popov determinant is x\~ m ~ x and the partition function becomes 



Z = J dx^-^e-^ 



(2) 



Due to the cancellation between fermions and bosons the partition function (2) coin- 
cides with the one that one would obtain from a purely bosonic n — m dimensional vector 
model. 

Let us consider now a supermatrix model defined by the partition function 



„ n+rr 

=/ n 

** A D 



n+m 

dA A B e~ s ^ (3) 

A,B=1 



where A is a supermatrix and it has the following block structure: 

{A- B } = whe re -■»= 1 - 2 - (4) 

l 4' b JV<J J a,P= l,2,---,n ^ ' 

where M\ andiV^g are even Grassmann variables while ip a p and ip b are Grassmann odd. In 
what follows we shall consider the case where A is a hermitian supermatrix and the action 
S[A] is invariant under the infinitesimal U(m\n) transformations 

5 A = [e, A] (5) 

where e is an antihermitian supermatrix and hence a generator of the supergroup U(m\n). 
The action S can be a function only of the U(m\n) invariants: 

str(A k ) = X)("1) TA ( Ak ) A A = E( A ") a * " E( Afc )°« 

A a a (6) 

a a 

where A a and /U a are the eigenvalues of A and also form a complete set of independent 
U(m\n) invariants. 

The partition function Z however, as given by the r.h.s. of eq.(3), is not well defined 
as the integral either gives identically zero or is in fact divergent. This can be seen directly 
in the m = n = 1 case where A is given by 



and the eigenvalues of A are: 



A = A + 



(8) 



A — fx 



A — /j, 



After denning £ + = A + fx and £_ = A — fx the partition function can be written as: 



The integral at the right hand side of eq.(9) can diverge at £_ = or, if £ vanishes at 
£_ = 0, at infinity in the (£+, £_) plane. This is the case for instance of any action of the 
type 



On the other hand, if none of these divergences occur, Z vanishes as the integrand is a 
total derivative in £ + . This latter property is a consequence of supersymmetry and it is 
completely general: the integral at the r.h.s. of eq.(3) involves the integration over 2nm 
fermionic variables while S(A) is a function only of the eigenvalues of A. The divergence 
of the integral in eq.(3) for any action of the type given in eq.(10) is also a property valid 
for the general U{n\m) model; in this respect supermatrix models are very different from 
ordinary matrix models where the integral defining the partition function is well defined at 
least for actions analogue to (10) with even I and suitable coupling constants. Besides, in 
matrix models the perturbative expansion of Z in powers of (k > 2) can always, at least 
formally, be defined as its coefficients are given in terms of convergent integrals, while in 
supermatrix models such coefficients are divergent as already the quadratic term in (10) 
in not positive definite. 

In order to define Z in a meaningful way we first reduce the integral at the r.h.s. of 
(3) to an integral over gauge invariant quantities, i.e. the eigenvalues of A. In doing that it 
is essential to preserve the U(n\m) invariance , so we fix the gauge in which A is diagonal 
and compute the corresponding Faddeev- Popov super determinant. In the chosen gauge 
the variation of A under an infinitesimal supergauge transformation is given by: 



where Xa are the eigenvalues of A. The Faddeev-Popov superdeterminant is then given 




(9) 




(10) 



k=2 



5A A B = e A B (X B - X A ) 



(11) 



by: 



sdet 




2X(-1)^4 +T B 



lla>b(Aa-A b ) 2 n a >/3(^a-^) 2 



(12) 



where {A^} = {A a ,// a }. 
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Eq. (12) is the supermatrix generalization of the volume element found by Mehta for 
ordinary hermitian matrices [4], and it reduces to it for m = 0. It follows that after fixing 
the gauge the partition function (3) for hermitian supermatrices becomes: 

71 \ f TT A\ TT A n a >b( X a-\b) 2 Ylaypi^-f^p) 2 -S\X a ,a a ] /-i n\ 

Z(n,m) = j I d\ a I I d/i a = ^ e * l a ' MaJ (13) 

By taking the volume element to the exponent one sees that the partition function (13) 
describes a one-dimensional Coulomb gas of opposite charges whose positions are given by 
the two sets of eigenvalues A a and \i a . The original action 5 acts as an external field. The 
integral in (13) has two kind of divergences: one, corresponding to the vanishing of the 
denominators in the integral, denotes the instability of the system towards annihilation of 
opposite charges, the other, which occurs for instance if 5 is given by (10), is the result 
of the instability originated by the external field that attracts the appropriate charges to 
±oo. The first type of divergences can be regularized by restricting the integration volume 
in eq. (13) to the region defined by the inequalities |A — /i Q | < e. In the limit e — > the 
integral is dominated by the configurations where the maximum number of denominators 
vanish. More precisely we find : 

Z(n,m) oce- m Z(n-m,0) + O(e- m+1 ) (14) 

where we assumed n > m and Z(n — m, 0) coincides with the partition function of a 
model of ordinary (n — m) x (n — m) matrices. Consequently, the correlation functions 
of the U(n,m) invariant operators O r = strA r coincide , in the e — > limit, with the 
corresponding correlation functions of an ordinary matrix model. The divergences that 
arise from the potential not being bounded can be eliminated by replacing for instance the 
action (10) with 

l n m n m 

S(X a , n a ) = 9k [E A* - e ^ ] + v [E x « p + E ^} ( 15 ) 

k=2 a=l a = l a=l ct=l 

where r\ > and 2p > I. In the limit rj — > and , for instance, / even the integral is 
dominated by the configurations where one set of variables, let's say the A's, are very large 
in absolute value (order rj~^p). 

Clearly the two limits e — > and rj — > do not commute, and the order of the two 
limits has to be given to define the theory. However the theory defined by first taking the 
limit e — > is the one that closely resembles ordinary matrix models. In fact if one performs 
a perturbative expansion of the partition function (3) in terms of Feynman diagrams, as 
in ref. [2] , the only divergent quantity one encounters is the propagator. This needs to be 
regularized by introducing a cutoff rj as in (15), or by any other equivalent regularization 
procedure (in ref. [2] for instance a Wick rotation in the /it's is performed ). Such cutoff 
is kept to the end of the calculation. On the other hand by performing the integral (3) 
over the matrix elements without fixing the gauge one takes into account automatically 
that configurations where couples of A and \x eigenvalues coincide have an infinite weight 
compared to the other configurations, and the e — > limit is automatically done. 
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Our analysis confirms the results of ref. [2] that supermatrix models are equivalent 
to the ones based on ordinary matrices; however this conclusion should be restricted to 
models where the action, as in eq (10), is given by a linear combination of terms of the form 
strA k . For more complicated potentials supermatrix models may differ substantially from 
ordinary matrix models. The supersymmetric extension of Penner-type models [5] could 
be an interesting example of this situation. Consider for example a supermatrix model 
coupled with r complex supervectors .As a result of the integration on the supervectors 
one gets in the integral at the r.h.s. of eq.(3) an extra factor sdet(A r ) which enhances the 
configurations with large A's and small /u's . This effect can break the pairing of eigenvalues 
due to supersymmetry and prevent a complete cancellation between bosonic and fermionic 
excitations. 

As a final remark it is worth noticing that, as in ordinary matrix models, different 
choices of A lead to different models:for instance if A has the form: 

A^ = g AC A CB (16) 

with 

Acs = {-1) TCTb Abc (17) 

and 

K S } = (C A) < 18 ) 

with C al3 a symplectic metric, then the action is invariant under transformation (5) with 
e belonging to the algebra of the orthosymplectic group and the corresponding Faddeev- 
Popov superdeterminant is given by: 

n a >jA a -A b | n a>/ 3(fa-^) 4 ( s 

EMAa-^a) 2 1 ) 

In the last equation the greek indices run from 1 to m/2 to take into account the fact that 
each eigenvalue \i a stands for two coincident eigenvalues of the symplectic part of A. 
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